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The smooth evolution of the tunneling gap of 
Bi2Sr2CaCu2 08 with doping from a pseudogap state in the 
underdoped cuprates to a superconducting state at optimal 
and overdoping, has been interpreted as evidence that the 
pseudogap must be due to precursor pairing. We suggest an 
alternative explanation, that the smoothness reflects a hidden 
SO(N) symmetry near the (tt,0) points of the Brillouin zone 
(with N = 3, 4, 5, or 6). Because of this symmetry, the pseu- 
dogap could actually be due to any of a number of nesting 
instabilities, including charge or spin density waves or more 
exotic phases. 

We present a detailed analysis of this competition for one 
particular model: the pinned Balseiro-Falicov model of com- 
peting charge density wave and (s-wave) superconductivity. 
We show that most of the anomalous features of both tun- 
neling and photoemission follow naturally from the model, 
including the smooth crossover, the general shape of the pseu- 
dogap phase diagram, the shrinking Fermi surface of the pseu- 
dogap phase, and the asymmetry of the tunneling gap away 
from optimal doping. Below T c , the sharp peak at Ai and 
the dip seen in the tunneling and photoemission near 2Ai 
cannot be described in detail by this model, but we suggest 
a simple generalization to account for inhomogeneity, which 
does provide an adequate description. 

We show that it should be possible, with a combination 
of photoemission and tunneling, to demonstrate the extent 
of pinning of the Fermi level to the Van Hove singularity. 
A preliminary analysis of the data suggests pinning in the 
underdoped, but not in the overdoped regime. 



I. INTRODUCTION 

A. Precursor Pairing? 

Recent photoemission0~i and tunnelingJ3~t] studies 
have provided a picture of unparalleled detail of the open- 
ing of the pseudogap in the underdoped cuprates. The 
most remarkable feature is that the pseudogap evolves 
smoothly into the superconducting gap as doping in- 
creases. This has led a number of researchers to conclude 
that the pseudogap must itself be related to superconduc- 
tivity: that it represents a form of short-range supercon- 
ducting order, or precursor pairing. It is the purpose 
of the present paper to show that this is not a foregone 
conclusion: there is an alternative interpretation (better: 
class of interpretations) in which the pseudogap repre- 
sents a competing ordered state. In this case, the appar- 
ently smooth evolution is due to an underlying symmetry 



of the instabilities of the nmblem -..a, manifestation of an 
SO(N) group, with N = M, 4B, SB, or ail 

For such an interpretation to hold, certain very strict 
conditions must be met. Specifically, there must be a Van 
Hove nestingo, with the Van Hove singularity (VHS) 
pinned close to the Fermi level ep over an extended dop- 
ing range. This prediction is now within the realm of 
experimental test, and a preliminary analysis of the ex- 
isting data seems to confirm the pinning. 

To describe the competing order parameters, we ana- 
lyze p_simple generic model, the pinned Balseiro-Falicov 
(BF)LLa model. Within this model, the total gap has a 
maximum at (jr, 0) in the Brillouin zone, given by 



r>2 



(1) 



where is an (s-wave) superconducting gap and Gk a 
charge-density wave (CDW) gap, defined below. This is 
exactiK the form proposed phenomenologically by Loram, 
et al.ta, and it immediately explains the smooth evolu- 
tion of the gap with doping: there is a single gap (at 
(7r, 0)) at all dopings, even though the system changes 
over from CDW near half filling to superconductor near 
optimal doping. The form of Eq.|l| immediately follows 
from an approximate SO (4) symmetry. It is in principle 
possible to disentangle the nature of the gaps from their 
behavior away from (tt, 0), although this involves sym- 
metry breaking terms, and hence is considerably more 
model dependent. The SO(N) symmetry also suggests 
that the pinned BF (pBF) model should provide a good 
approximate representation for a wide range of compet- 
ing phases - in particular the pseudogap phase may also 
be an antiferromagnetic phase or a flux phase (which is a 
form of dynamic CDW) , or indeed a striped phase which 
is a combination of two of these phases. 



B. Van Hove Pinning 

In the generalized Van Hove scenarioEll, there are two 
separate phenomena which contribute to Van Hove pin- 
ning. First, as part of the Mott-Hubbard transition, 
strong Hubbard-U correlation effects renormalize Cu-0 
hopping t to zero at half filling, leaving a residual energy 
dispersion associated with exchange J. In the absence 
of correlation effects, the 'bare' VHS would fall at a fi- 
nite hole doping fixed by the band parameters (t' in the 
tt'J model) - this doping will be close to, but not nec- 
essarily the same as optimal doping. However, since sec- 
ond (Cu) neighbor exchange is expected to be small, the 
exchange bands have the simple dispersion J(c x + c y ), 
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where Cj = cosfcia. Hence, at half rilling, the Fermi level 
approximately coincides with the VHS. This phase can 
be further stabilized by Van Hove nesting, which opens 
up a large gap in the dispersion near (tt, 0). A good fit. 
to the dispersion in the related insulating compoundtS 
Sr 2 Cu0 2 Cl 2 (SCOC) can be founcM~El by -assuming that 
the nesting is associated with a flux phaseEa. 

With doping, t is gradually restored, producing a be- 
havior which cannot be described by a rigid band filling 
model. At first, the VHS shifts faster than the Fermi 
level, so the VHS lies below the Fermi level, but very close 
to it. Gradually, t saturates, the VHS stops shifting (usu- 
ally close to its bare value) , and at some finite doping x c 
the Fermi level crosses this new VHS. Thus, the Fermi 
level coincides with the VHS twice: at half filling and 
near the bare VHS, and remains anomalously_close at in- 
termediate dopings. Since its initial discoverycij, this cor- 
relation induced, pinning has been confirmed by a number 
of calculations^. 

But even stronger pinning is possible. Near x c , the 
energy can be significantly lowered by a second nesting 
instability. Assuming this second instability to be charge- 
density wave (CDW) related, a self-consistent 3-band 
slave boson calculation demonstrated that this model can 
lead to two free-energy minima, one at half filling and 
the other at x c Eil. This results in striped phases, with 
each phase pinned near a VHS: the magnetic stripes near 
the J-dominated VHS near half filling, and the charged 
stripes near the t-dominated VHS at x c . Since the stripes 
are nanoscale (due to long-ranged Coulomb repulsion), 
the system evolves rather smoothly with doping, with 
the Fermi level remaining even more strongly pinned to 
the VHS at all dopings. 

While the above model is in good agreement with 
experiments on both the pseudogaps and the striped 
phasesE3, the identification of the specific nesting insta- 
bility phases is less secure. This is because the nesting 
and pairing instabilities of the Van Hove scenario form 
a group, SO (6), and the instabilities associated with twa 
dual 6-dimensional superspins are nearly degeneratetJ, 
so the issue of which phase is the most unstably de- 
pends sensitively upon secondary parameters which are 
not well known. The various possibilities include the an- 
tiferromagnetic and d-wave superconducting instabilities 
of Zhang's SO(5)E3, as well as CEjW's, flux phases, and 
a more exotic spin current phaseES. (Two-leg ladders 
havei-an even larger assortment of instabilities to choose 
frorrO.) It is difficult to incorporate the details of this 
phase separated regime rigorously into calculations of the 
ARPES and tunneling spectra. 

Nevertheless, a remarkably simple picture of the pseu- 
dogap (VH nesting) phase diagram can be developed via 
a simple Ansatz for the pinned striped phases. A one 
band model is assumed 

e fe = -2t(c x + Cy) - 4t'c x c y (2) 

(with Ci = cosfcia), and the second neighbor hopping is 
adjusted with doping to pin the VHS to the Fermi level 



over a doping range, from half filling x — to a doping 
x c . ApproximatelycH, 

2t' 

t = — = - 1.04 tanh 2.4a;. (3) 

Within this model the striped phase is represented by a 
single nesting instability which splits the band dispersion 
at (7r, 0). Such a model was initially introduced (with- 
out the pinning) by Balseiro and Falicovll3 to study the 
competition between CDWs and s-wave superconduc- 
tivity. We have employed this pinned Balsepp-Falicov 
(pBF) model in previous pseudogap studiesE3~E2l and will 
continue to use it here. [It must be stressed that we do 
not employ the CDWs to mimic the spatial pattern of 
the stripes, but rather to reproduce the Fermi level pin- 
ning.] We are currently generalizing the model to include 
d-wave superconductivity and a variety of other nesting 
instabilities. In the present version, there are electron- 
phonon coupling energies \q associated with CDWs, and 
Aa with superconductivity, which may or may not be 
equal. The gap equations are solved self-consistently, 
with free parameters t, x c , the A^'s, and u) p h, a bosonic 
cutoff frequency. There is relatively little data in the 
overdoped regime, so we choose a simple model: for 
x > x c , the band parameters cease evolving, and the 
additional holes just rigidly fill the band, shifting the 
Fermi level away from the VHS. This simple picture ap- 
pears capable of describing the overdoped regime in both 
YBa 2 Cu 3 07_<5 (YBCO)Tvhere the pseudogap vanishes 
close to optimal dopingo and La 2 _ K Sr a ;Cu04 (LSCO), 
where the pseudogap appears to persist well into the over- 
doped regime^j. However, a more complicated behavior, 
including a second range of two-phase coexistence (Sec- 
tions 11.6-11.7 of Refill!) , is not ruled out. 

Within this picture, there is a natural doping depen- 
dence associated with the competition between nesting 
and pairing instabilities. At half filling, t' — 0, the per- 
fect nesting overwhelms the pairing instability, leading to 
a pure nesting instability. As t' increases with doping (to 
maintain the VH pinning), the nesting gets worse, while 
pairing is less affected (it is actually enhanced, Section 
II. B). This leads to a crossover to a pairing instability as 
a function of doping. Since nesting does not gap the full 
Fermi surface (when t' is large enough) , superconductiv- 
ity can appear at a lower temperature, with a T c which in- 
creases with doping. As long as the VHS remains pinned 
to the Fermi level and the electron-phonon coupling re- 
mains doping independent, the nesting T n will decrease 
with doping x while the superconducting T c will increase. 
At some point, the two transitions would cross. However, 
at this point T n will be rapidly suppressed to zero, since 
superconductivity gaps essentially the full Fermi surface 
(the small residual Fermi surface for a d-wave gap could 
only sustain a nesting instability at a much lower T). 
While exactly this behavior is found in YBCO, in LSCO 
and probably Bi 2 Sr 2 CaCu 2 0s (Bi-2212) as well the pseu- 
dogap persists into the overdoped regime. This behavior 
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can be modelled by having the strength of the bosonic 
pairing decrease with increasing x. 

Since nesting splits the VHS degeneracy, we introduce 
some notation to clarify the following discussions. The 
lower VHS (VHL) is the VHS shifted below the Fermi 
level, at energy El, and hence visible in photoemission. 
The upper VHS ( VHU) is shifted to Eu , above the Fermi 
level, and hence can only be seen in tunneling (or inverse 
photoemission). In conventional plots of tunneling spec- 
tra, VHL appears at negative voltages, in the electron 
extraction mode, and VHU at positive voltages. The Van 
Hove centroid (VHC) is the average position of these two 
features: Ec = (Eu +El)/2 - it is the position at which 
the single, unsplit VHS would fall at high temperatures, 
well above the pseudogap transition. 

In terms of these features, the photoemission gap is 
Ap e — £f — El- We will show (Fig. || below) that the 
tunneling gap is Atu — {Eu — El)/2. In this case, 
the pinning of the Fermi level at the VHC, ep ~ Ec, 
can be rewritten in terms of measurable quantities as 
A TU = A PE . In Section III.C, below, we will demon- 
strate that this relationship appears to be satisfied in the 
underdoped regime. 



II. STRIPED PHASES AND PSEUDOGAP 

A. Phase Diagram of Pinned BF Model 

For completeness, we recall the energy disper- 
sion and the gap equations of the BF model. In 



terms of the four-component wave vector 

i^p ^4+<5> ^-fc' "^-(fc+Q))' the mean neld BF Hamil 
tonian tensor is 



Hbf = 
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In terms of a function 



1, if | eg - e F \ < hLO ph ; 
0, otherwise, 



(4) 



(5) 



the gap functions are = A0£ for superconductivity, 
and G £ = G + Gi&^Q^ for the CDW. The energy 
eigenvalues are -E±,fc and their negatives, with 
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l+Q) 

e% + + (A fe - A k+Q ) 2 , e k ± = e k ± £/c+q, and the nesting 
vector Q = (n,Tr). If the magnitudes of the (attractive) 
phonon-induced electron-electron interaction energies are 
Aa and Ag, then the gap equations are 
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with Go = ©fc©fe + Qi ©i — 1- A similar modeled has 
recently been applied to analyze the photoemission asso- 
ciated with a pure CDW phase. 

In the pinned BF (pBF) model, the Fermi level is 
pinned to the VHC, via Eq. ^, for doping x between 
half filling (x — 0) and some critical doping x c , while 
for x > x c , the curvature is fixed at t' = t'(x c ), while the 
Fermi level shifts off of the VHC in a rigid band filling. 
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FIG. 1. Phase diagram of pinned Balseiro-Falicov model. 
Circles = net tunneling gap, A t ; solid line = CDW transition 
temperature T p ; dashed line = superconducting transition T c ; 
dotted line = T c in absence of CDW, x c > 0.5. 

For fixed values of the parameters t, x c , Aj, and phonon 
cutoff uj p h , the pseudogap phase diagrams are derived 
by solving Eqs. self-consistently. For Fig. [|, the pa- 
rameter values are chosen as t = \\ = Xq = 0.25eU, 
w p h = 45meU, and x c — 0.123. At half filling, per- 
fect nesting wins out over superconductivity, but with 
doping the nesting becomes poorer, while superconduc- 
tivity is enhanced, leading to a crossover. As soon as 
the superconducting T c is larger than Tcdw, Tcdw is 
rapidly suppressed to zero. The dotted line in Fig. [j] 
shows how T c would evolve in the absence of the CDW, 
if the Fermi level remained pinned to the VHC over the 
full doping range (x c > 0.5). A typical temperature de- 
pendence of the resulting gaps is illustrated in Fig. ||, for 
x = 0.1. Note that A t , Eq. [l2|, is found self-consistently 
to evolve smoothly with temperature; this implies that 
Gfc must actually decrease when superconductivity ap- 
pears, A k ^ 0. 
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FIG. 2. Evolution of gaps in pinned Balseiro-Falicov model. 
Solid line = superconducting gap A k ; dashed line = CDW gap 
Gk = Go + Gi; dotted line = net gap, A t , Eq. Parameters 
correspond to Fig [l|, x = 0.1. 

To reproduce the phase diagram of YBCO, the 
crossover must arise close to Xc, as in Fig. |, while m 
LSCO it falls well after x c , Fig. y. The phase diagram of 
LSCO can be modelled by choosing Ac 7^ Aa, and letting 
the latter vary with doping, Fig. H. 
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FIG. 3. Model pseudogap phase diagram for LSCO. Pa- 
rameters x c — 0.16, Xq = O.GeV, with doping dependent Aa 
(inset). Solid line = CDW transition temperature Td\ dashed 
line = superconducting transition T c . 

One feature of Fig. || should be noted. In LSCO the 
pseudogap appears in-some experiments to persist into 
the overdoped regimeu'L-il, while other experimentp-fixid 
a striking crossover in properties at-aptimal dopingoLO. 
In the thermodynamic experiments^, the pseudogap ap- 
pears to be closing linearly with doping, but then has 
a break in slope at optimal doping, and falls off more 
slowly at larger doping. Figure mimics this behav- 
ior: the Fermi level is pinned at the VHC until optimal 
doping x c = 0.16, and then depins. Despite the fact 
that the Fermi level now shifts with doping away from 
the VHC, the CDW instability is strong enough that the 



CDW phase persists out to x = 0.27, at which point 
it becomes unstable (discontinuously). Such behavior is 
only possible because the superconducting transition is 
suppressed far below T p . 

The pBF model is a 2D mean field theory. In a more 
accurate three-dimensional calculations, the transition 
temperatures in the above phase diagram are likely to be 
replaced by crossover temperatures at which 2D fluctua- 
tions become strong - i.e., temperatures at which a pseu- 
dogap opens. When interlayer coupling becomes strong 
enough, a real three-dimensional order, with true energy 
gaps, can develop. 

A (weak-coupling) limitation of the BF model is the 
discontinuous change in A and G\ when the energy 
crosses uj p h- When this feature becomes prominent in 
the dispersion, the BF results must be considered as only 
a qualitative indication of the results of a proper strong 
coupling calculation. 



B. SO(6) 

Whereas the above discussion has been in terms of a 
particular competition - between a CDW and an s-wave 
superconductor, the results are indicative of a much more 
general situation. That is because the VHS. instability 
has an underlying SO (6) symmetry groupES. There is 
actually a pair of 6-dimensional 'superspins', which con- 
sist of various nesting or pairing instabilities of the VHS. 
Nesting (pairing) operators are those instability opera- 
tors which do (do not) commute with the number oper- 
ator Q. _ 

For example, Zhang's SO(5) grouprJ is a subgroup of 
this SO (6), where the competition is between antiferro- 
magnetism (or a spin-density wave - SDW) and d-wave 
superconductivity. The SO (6) offers greater flexibility, 
with the nesting instability being either the CDW-, or 
SDW, or a flux phaseca or a spin current instabilityc3. 
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FIG. 4. Susceptibilities \Q (dotdashed line) and \o (solid 
line) vs. band filling 1 + x for Eq. ^ with At' = Ef- 

The generic features of the instabilities can be most 
clearly seen be looking at the doping dependence of 
the corresponding (bare) susceptibilities yoXa, w). For 
a pinned VHC, this is illustrated in Fig. For all 

the instabilities, there are only two bare susceptibili- 
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ties, xo = Xo(<? = 0,0) for the pairing instabilities and 
Xq = Xo(<l = Q7O) for the nesting instabilities. At 
half filling these susceptibilities are degenerate, but as 
t' increases, the VHS's become more one-dimensional as 
the plaaes become more like interpenetrating Cu-O-Cu 
chainsEa. In this case, nesting is greatly reduced, since 
the chains run at right angles, whereas pairing is en- 
hanced, because the ID bands have a stronger VHS. 

We note that the phase diagram of Fig. [l] shows con- 
siderable evidence for the underlying SO(N) nature of 
the VHS: (a) despite a crossover from CDW to super- 
conducting, the evolution of the total gap A t and tran- 
sition temperature T t with doping is extremely smooth; 
indeed, the gap ratio 2A(0)/fcsT c ~ 4.1 is nearly in- 
dependent of doping, (b) The total gap is given by a 
vector addition, Eq. |l|. And (c) the overall shapes of 
the individual (decoupled) transitions Tj (x) resemble the 
doping dependences of the generic susceptibilities, Fig. ^. 
Thus, the solid line in Fig. p] (pseudogap phase) closely re- 
sembles the nesting susceptibility (a larger susceptibility 
corresponding to a higher transition temperature) , while 
the pairing susceptibility resembles the superconducting 
phase (dotted lines in Fig. |l|). 

Note finally that just these features - in particular fea- 
ture (a) - have been taken as evidence that the pseudogap 
must itself be related to superconductivity as a form of 
precursor pairing. Instead, we find that the phase dia- 
gram of Fig. U should look fairly similar for any compe- 
tition between a nesting and a pairing instability - e.g., 
either CDW vs s-wave superconductivity or SDW vs d- 
wave superconductivity. The actual phases observed will 
depend on the interaction parameters, and may vary with 
doping. Hence, it seems unlikely that the experimental 
pseudogap represents superconducting fluctuations, un- 
less all nesting instabilities are somehow strongly sup- 
pressed. 

We stress that the similarity between Figs. [I] and |] 
exists even though the BF model does not possess SO (4) 
symmetry. The form of the phonon coupling leads to very 
different gap equations for CDW and superconductivity, 
Eqs. H"^- Thus, the two instabilities are not degenerate 
at zero doping. Yet the doping dependence of the two 
instabilities follows the susceptibilities of Fig. [|, and the 
gap at (tt, 0) is just the total 'length' of the two individual 
gaps (adding in quadrature). _ 

Experimentally, there is evidence for striped phasesLJ 
at intermediate dopings. In SO(5), these stripes can be 
interpreted as a combination of magnetic and supercon- 
ducting stripes. However, in this case it is not clear how 
to interpret the 1/8 anomaly (x = 1/8) where the stripes 
show long range order but there is no superconductiv- 
ity. It seems more likely that superconductivity competes 
with the stripes, and that the stripes arise from a com- 
petition between two nesting instabilities, the flux phase 
at half filling and a CDW near xjE3. It is important to 
recognize that the flux phase does not involve the elec- 
tron spins - the magnetic moments are orbital. Hence, 



the phase is beat understood as a dynamic CDW phase 
(compare Ref.ufj), so a model which approximates the 
pseudogap by a CDW Ansatz (such as the pBF model) 
should provide a quite reasonable first approximation. 

C. Tunneling Spectra 

Under certain special conditions, both the tunneling 
and photoemission can reveal very direct information 
about the spectral function of the interacting elections. 
The quasiparticle tunneling current can be writtcreil 

/OC T 
,±A R (k,e)A L (p,e + eV)x 

[n F {e)-n F (e + eV)], (9) 

where is the tunneling matrix element, Ai is the ap- 
propriate spectral function in the metal on the left (L) 
and right (i?) of the tunneling junction. If the tunneling 
matrix element is considered to be constant, independent 
of k and p, it can be taken out of the integral, yielding 

r°° de 

I = 2eT 2 / —N R (e)N L {e + eV)[n F (e) - n F {e + el/)], 

(10) 

with Ni the appropriate dos. For an NIS junction, taking 
Nl to be energy independent, the tunneling conductance 
is 

dl e 2 T 2 

G=^- = N L N R (-eV). (11) 

The conventional wisdomB is that tunneling is not sen- 
sitive to the dos: the explicit factor of dos in Eq. |^ is 
cancelled by T 2 oc v Fl the Fermi velocity. However, Wei, 
et al.EJ have shown that this cancellation breaks down in 
the presence of strong anisotropy. In particular, tunnel- 
ing along the c-direction into a two-dimensional metal 
directly measures the in-plane dos (at least in the thin 
junction limit). Hence, in accord with Eq. [ll], the tun- 
neling conductance is proportional to the tunneling dos. 
This result is the basis for the present analysis. For three- 
dimensional materials, there have been previous propos- 
als that the VHS's should show up in tunnclingCj, and 
the tunneling spectra of CDW superconductors have been 
analyzedEJ. 

The tunneling and photoemission are derived from the 
spectral function of the model 

A(k,u) = 2nJ: i=± [ul k S(u;-E hk ) + vl k S(uj + E^ k )}, 

(12) 

with eigenenergies given by Eq. ^, and coherence factors 
u±,k = u±, o cos(0 + (j> T ), (13) 
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2 1 r^i.A; + £ fe £ fc+Q ^ A k&k+Q +G k 
COS 0± = -(1 + [ ■ ~— ]), 



sin 20 : 



sin 2<p± = 



(16) 



(17) 



(18) 



In deriving the spectral function, it is convenient to use 
the group, theoretical techniques of spectrum-generating 
algebrasE^I. 

For a pure CDW (Afc = 0), the spectral function sim- 
plifies to 

A(k, u) = 2n[u 2 k 6(u -E k+ )+ v 2 k S{uj - (19) 
with ul = l-v 2 k = (l + e k -/E' k )/2, 

E k± = (e fc+ ± E' k )/2, (20) 



andE' k = Jel_+AGl 



> 




r x s r n(e) 

FIG. 5. Energy dispersion for data of Fig. ^, with r = —0.2, 
T = 10K. Right = tunneling dos. (a) = full dispersion; (b) 
= blow-up of region near ef. 



Figure [| shows the energy dispersion for a com- 
bined CDW-superconductor (left) and the associated dos 
(right). Part (a) shows the full dispersion, while (b) is 
a blow-up of the region near the Fermi level. There are 
four bands; the CDW order folds the M-point of the Bril- 
louin zone into the T point, while superconductivity folds 
both of these bands around the Fermi level (dotted lines). 
However, these ghost bands carry little spectral weight 
(coherence factor << 1) except very near the Fermi level 
- note that there is no indication of the tops or bottoms 
of the dotted lines in the dos. Figure ^ replots the dis- 
persion, giving an indication of the spectral weight. See 
further Fig. |Io|, below. 




FIG. 6. Replot of energy dispersion of Fig. || illustrating 
spectral weight. Coherence factor > 0.6: solid lines; between 
0.1 and 0.6: dashed lines; between 0.001 and 0.1: dotted lines. 

A close look at the region near the Fermi level, 
Fig. ||b, reveals that structure in the tunneling dos is di- 
rectly related to features in the dispersion of the gapped 
bands. Thus, peak A is associated with the dispersion 
at (tt, 0) - the VHS peak split by the combined CDW- 
superconducting gap. Peak B is due to the superconduct- 
ing gap away from (tt, 0) - particularly near (7r/2,7r/2). 
Accordingly, it will be considerably less prominent for 
a d-wave superconductor, where the gap vanishes at 
(7r/2,7r/2). Feature C is associated with_the CDW gap 
Gk near (7r/2,7r/2). As discussed in RcfO, Equation |20| 
can be rewritten as 



E k + = -At'c 



± 



(21) 



showing that the CDW gap is fixed to the VHS, but not 
to the Fermi level, so that at fc-points away from (tt, 0) 
there can be two gaps aLxhfferent energies. Feature C is 
further discussed in RefO. Whether this is only a weak- 
coupling effect that would be washed out in a strong- 
coupling theory remains to be seen. Finally, feature D is 
associated with the phonon-related discontinuities in the 
dispersion at uJ v h- 
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The appearence of two gap-like features in the tun- 
neling spectrum is not necessarily a consequence of two 
competing order parameters, but of gap anisotropy. In- 
deed, for a pure (generalized) s-wave superconductor, 
there will be two peaks in the tunneling dos whenever 
the gap is anisotropic, as long as the minimum gap is 
non-zero. Moreover, even for a pure d-wave supercon- 
ductor, there can be two gap-like features if the VHS is 
not at the Fermi level, Fig.[lO|. 

D. Line Broadening 

In the one-dimensional CDW, the principal source of 
line broadeping-, above the Peierls transition is CDW 
fluctuationsoeS. The phonon propagator diverges at 
the transition: 

D(q) * * (22) 

i~ 2 + W-Q\ 2 

with oc (T — T p ). For the tap-dimensional Van Hove 
problem, a similar form holdsE3 with T p = in the ab- 
sence of interlayer coupling. This results in an intrinsic 
broadening T — vp oc y/T. 

In the superconductor, a true long-range ordered state 
is possible, with a real gap at the Fermi surface leading 
to a considerably reduced T. Such a reduced scattering 
has been observed in a number of transport properties, 
as well as in the photoemission spectrum. However, the 
scattering should be restored at high frequencies, when 
the CDW fluctuations can break a rpalr. On analogy 
with the results of Coffey and CoffeyCJ, we assume that 
the crossover will fall near 2A&. (Strictly speaking, the 
crossover should be at 3A^ for an s-wave superconduc- 
tor.) 

III. COMPARISON TO EXPERIMENT 

A. Experimental Situation 

Both photoemission and tunneling measurements of 
the pseudogap have been largely restricted to Bi-2212. 
While the general features of the spectra are reasonably 
well understood, there are still differences in detail. In 
this subsection, we summarize a number of unresolved 
issues. 

(a) In optimally-doped and lightly underdoped Bi- 
2212, photoemission detects a significant rearrangement 
of spectral weight in the superconducting state. The low 
temperature (T < T c ) spectra are generally characterized 
by three features: a sharp quasiparticle peak Ai sepa- 
rated from the Fermi level by a small (~ 25-50 meV) 
gap and a broad hump at much higher energies, A2 sep- 
arated by a well-defined dip at energy ~ 2Ai. Above T c , 
the peak and dip disappear, leaving a feature similar to 
the hump. 



The dip feature is now widely interpreted as coupling 
to some form of collective mode, the exact nature of 
which is not well understood. The fact that it varies with 
doping, scaling as 2Ai is puzzling. One possible interpre- 
tation is that it is related to electron-electron scattering. 
The broadening of the photoemission features is known 
to decrease dramatically just below T c , and this decrease 
would be expected to fall off at frequencies much larger 
than the gap (at 3A for an s-wave superconductor, 2A 
for d-wave). Such a result could be very useful in sepa- 
rating the superconducting gap from the pseudogap, but 
fits to this model have been unsuccessful in reproducing 
the dip amplitude. We will show (Section III.E) that it is 
possible to get a large dip if there is significant spectral 
weight in the hump feature. 

(b) The presence of the dip has led to some confu- 
sion on the doping dependence of 'the gap': should one 
count Ai ox A2. The latter choice was made by Mar- 
shall, et al.p. In their earlier work, the Argonne group 
used neither choice, but defined the gap in terms of the 
leading edge near ep. This feature should scale with the 
peak position of Ai, but be somewhat smaller. In their 
more recent work they have used the peak position Ai for 
better comparison with tunneling studies. Our analysis 
suggests that this is a more appropriate choice, and we 
shall follow this latter usage. Our own belief is that, since 
the hump persists above T c , A2 is a proper candidate for 
the pseudogap. However, it is not always easy to extract 
this feature from the published literature. Fortunately, 
Dingo is now carrying out a careful study of both gaps. 

Near optimal doping, Ai ~ A2, so the literature un- 
certainty will not greatly affect our reconstructed phase 
diagram, Fig. [j], below. For stronger underdoping, the 
sharp peak and dip gradually wash out, while the hump 
shifts to higher energy. We find that this shift is clearly 
revealed in the doping dependence of T*, the tempera- 
ture onset of the pseudogapEl, and that these data are 
consistent with the Stanford measurements of the hump 
data, suggesting that the photoemission provides a single 
pseudogap phase diagram. Moreover, this diagram is re- 
markably consistent with the pseudogap phase diagrams 
found for LSCO and YBCO on the basis of transport 
measurements . 

(c) Tunneling is mainly sensitive to the sharp peak and 
dip feature in the superconducting state, and it is not al- 
ways clear whether the hump feature is seen. However, 
some groups clearly do see structure above T c El. One 
problem is that tunneling is very surface sensitive, and 
the BiO layer on the top surface has an insulating gap. 
However, consistency with the photoemission would sug- 
gest that a pseudogap. should be present above T c . Re- 
cently, Renner, et al.o have reported that the tunneling 
gap inside a vortex core resembles that found in the pseu- 
dogap above T c . 

(d) Thexe is a question about overdoping, with Ren- 
ner, et al.tj detecting the pseudogap in overdoped sam- 
ples, while photoemission studiesEj find a rapid collapse 
of the pseudogap in overdoped Bi-2212. It seems likely, 
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however, that in LSCD. the pseudogap persists well into 
the overdoped regimeEl 

(e) There is also some disagreement on the gap ratio 
2A(0)/fc B T*, with Oda, et a\li reporting a value 4-5, 
whereas a value ~ 8 can be extracted from the data of 
Ding, et al.B. 



agreement with the pseudogap determined from trans- 
port in LSCO and YBCOtii, Fig. §. In all three cases, 
the doping axis had to be scaled to cause the curves to 
coincide. The scaling suggests that, if the optimal doping 
for superconductivity in LSCO is x — 0.16, it is an effec- 
tive x = 0.2 in YBCO, and x = 0.32 in Bi-2212. Such a 
shift in x c is consistent with the Uemura plotH. 



B. Pseudogap Phase Diagram 

Figure |?] summarizes the photoemission and tunneling 
data on the pseudogap in Bi-2212. Most of the nhotoe- 
mission data are taken from the Argonne grouptj, but 
the two lowest dopings are from the Stanford, groupEJ, 
and from the insulating phase, S^CuC^CbtS. Since 
there has been disagreement on how to determine the 
pseudogap peak position, it is important to note that 
the Argonne data are consistent with a large increase 
in T* (open circles) at low doping. Figure g illustrates 
an optimized fit of the pBF model to the Bi-2212 pseu- 
dogap phase diagram. Details about choosing the pa- 
rameters are discussed in Appendix A. Parameters are 
t = 0.25eV, x c - 0.335, uj ph = 35m e y, A G = U6meV, 
A A = IMmeV. 




x 



FIG. 7. Pseudogap phase diagram in Bi-2212 determined 
from phetpemission and tunneling. Dipmonds = photoemis- 
sion gaptj; + 's = photoep«ission gapuiij; circles = mea- 
sured from photoemissionEa; squares = tunneling gapO. Solid 
line = guide to the eye. 

The overall qualitative and quantitative agreement is 
quite good. A number of features of the experimental 
data are worth commenting on. First, there is an approx- 
imately constant ratio between the total pseudogap, de- 
fined as the energy shift between the dispersion at (tt, 0) 
and the Fermi level, and the pseudogap onset tempera- 
ture T*, 2A(0)/fc s T* ~ 8. A similar but smaller ratio, 
4.1, is found in the calculation inote that this ratio is 
close to that found by Oda, et al.H). Secondly, the overall 
shape of the curve, strongly ^-dependent in the under- 
doped regime, nearly constant in the overdoped regime, 
is well reproduced by the theory. Thirdly, the pseudogap 
determined from photoemission in Bi-2212 is in excellent 
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FIG. 8. Model pseudogap phase diagram for Bi-2212. Solid 
line = CDW transition T p ; dashed line = superconducting 
transition T c ; circles = total gap At at IK. 
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FIG. 9. Pseudogap phase diagram in LSCO and YBCO 
determined from transport. Open circles: from resistivity of 
YBCO; other symbols: transport measurements in LSCO, see 
Ref£j; solid line from Fig.^, with x axis scaled by a factor of 
2. 

Finally, as noted previously!, there is excellent agree- 
ment between the photoemission pseudogap at (jr, 0) , 
Ape, and the tunneling pseudogap, Atu, defined as half 
the splitting between the two tunneling peaks. As illus- 
trated in Fig. |j| this result is also in excellent agreement 
with theory. Note that the theory holds specifically when 
the Fermi level is pinned at the VHC. Indeed, we believe 
that the present dataEl constitute the strongest proof for 
this pinning. To illustrate the strength of the evidence, in 
the following subsection we will test the null hypothesis: 
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assume that VH pinning is absent, and see how different 
Ape and Atu would be. 



C. Van Hove Pinning vs d-wave Superconductivity 

It has been suggested that what we have interpreted as 
Van Hove pinning can be alternatively explained as due 
to simple d-wave superconductivity, since the d-wave gap, 
Ad = Aq(c x — c y ) is largest near (tt, 0). Here, we demon- 
strate that this is not the case. We assume that dop- 
ing is accomplished by the filling of a rigid band, Eq. g, 
and there is a fixed d-wave gap with doping-independent 
strength 2Ao = 50meV. Figure [l(] illustrates the energy 
band dispersion near the Fermi level, and the correspond- 
ing tunneling spectra, at a series of dopings away from 
the VHS. The solid lines are the principal branches of 
the dispersion, with coherence factors > 0.5; the dashed 
lines are those parts of the ghost branches with coher- 
ence factors in the range 0.1-0.5. A large phonon energy 
Tiijph — 90meV was assumed, to shift the phonon peaks 
in the tunneling spectra out of the range of interest. 
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FIG. 10. Energy dispersion near the Fermi level (left), and 
associated tunneling density of states N(E) (right) for a pure 
d-wave superconductor, as the VHS sweeps through the Fermi 
level: (E V hs - E F )/t = -0.2 (a), -0.1 (b), (c), +0.1 (d). 

The tunneling spectra clearly show the VHS moving 
through the Fermi level. Note that the VHS is always at 
(tt, 0), but the superconducting gap is at the Fermi level, 
and hence in general away from (tt, 0). Superconductivity 
and the VHS actually produce two separate peaks in .the 
tunneling spectra. (This fact has been noted earlierE3.) 
If the spectra are too broadened to resolve the individ- 
ual peaks, one would find the lower (upper) peak to be 
more intense in underdoped (cwerdoped) samples. This 
behavior is seen experimentallyffl, and constitutes strong 
evidence that at optimal doping the Fermi level exactly 
coincides with the VHS - in agreement with an earlier 



predictionEZI. We claim additionally that the Fermi level 
shifts away from the VHS anomalously slowly on the un- 
derdoped side. 

It can be seen that the dominant gap arises away from 
the VHS - at the point where the bands cross the Fermi 
level. However, there is a subsidiary structure associ- 
ated with the saddle point dos peak, which rapidly shifts 
away from the Fermi level with doping. This would lead 
to a double peak structure in the tunneling spectrum. 
Moreover, it would mean that the tunneling gap (taken 
as the lower-energy, more prominent feature) is distinct 
from the photoemission gap (defined by the dispersion at 

(tt.0)). 
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FIG. 11. Doping dependence of Atu the tunneling (circles 
and dashed lines) and Ape, the (tt, 0) photoemission (x's and 
solid curve) gaps, for a d-wave superconductor in the absence 
of VH pinning (r = -0.25). 

This is illustrated in Fig. [n], which shows the doping 
dependence of these two gaps. The symbols were esti- 
mated from the curves of Fig. [l^ while the curves are 
analytical expressions, given by the curvature of the dis- 
persion relations (Appendix B). 




FIG. 12. Normalized splitting of Fermi level from the VHS, 
measured as the normalized difference between the tunneling 
and photoemission pseudogaps, (Ape — Atu) /Atu- Solid 
line = theory in the absence of pinning (similar to Fig. [H], 
but with t = —0.38); open circles = derived from data of 
RefsM 

In Figure we plot the calculated difference between 
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the tunneling and photoemission gaps, normalized to the 
tunneling gap, and compare this to our estimate of the 
experimental differences, Fig. ^. For convenience, we re- 
plot the tunneling and photoemission gaps from Fig. (?] in 
Fig. [l3| Only three data points could be extracted, and 
there are considerable error bars in the measurements. 
(Note in particular that Ape — A-tu cannot be negative.) 
The errors are smallest for the underdoped system: both 
tunneling and photoemission data are reported in Ref.H, 
and are presumably for similar samples. In the other two 
systems, tunneling data frorn— Ref.B are compared with 
photoemission data from Ref.til. 
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FIG. 13. Comparison of photoemission (diamonds)Ej and 
tunneling (squares)u gaps in Bi-2212. Solid and dashed lines 
= guides to the eye. 

Despite these limitations, the results are intriguing. 
There is a hint that the VHS depins from the Fermi level 
in the overdoped regime, but there is no sign of depin- 
ning in the underdoped samples. Thus, the data appear 
to rule out the hypothesis of rigid band filling, strongly 
implying that below optimal doping the Fermi level is 
pinned to the VHS. Clearly, this experiment needs to be 
repeated much more carefully, and over a wider doping 
range. In particular, (a) the photoemission gap must be 
measured exactly at (71", 0), and (b) both photoemission 
and tunneling must be carried out on the same (or iden- 
tical) samples, to minimize sample-to-sample variations. 
However, a direct determination that the Fermi level is 
pinned to the VHS would confirm a number of strong cor- 
relation theories, and would have a profound influence on 
future theoretical modeling. 

D. Tunneling Spectra 

Figure |lj illustrates the evolution of the calculated 
tunneling spectra with temperature in the underdoped 
regime (x = 0.274, dotted line in Fig. ||). The dotted lines 
show the normal state VHS above T p ; the dashed lines are 
in the CDW phase; and the solid lines are in the mixed 
CDW-superconducting phase. The only broadening is 
thermal, due to the Fermi function; the disorder broad- 
ening T = 0. As noted above (Fig. there are promi- 



nent features at the phonon frequency, uj p h = 3bmeV 
(arrows). For a realistic phonon spectrum, they would 
split up into the strong-coupling factor a 2 F. There have 
been a number of reports in the literature of the observa- 
tion of such structure, but it is not clear how reproducible 
it is. 
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FIG. 14. Tunneling dos for Bi-2212, using parameters of 
Fig.^ (dotted line). From bottom to top, the temperatures 
vary from 10K to 100K, in 10K intervals, then go to 120, 150, 
200, and 300K. Dotted lines, T > T v = 116-fT; dashed lines, 
T P >T >T C = 69K; solid lines, T < T c . All curves are offset 
for clarity (all essentially coincide for |eV| > lOOmeF). 
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FIG. 15. Tunneling dos for Bi-2212, using parameters of 
Fig.[t4| but with hu) p h=90meV . From bottom to top, the tem- 
peratures vary from 10K to 100K, in 10K intervals. 

In order to better see the shape of the tunneling peaks, 
the dos is recalculated in Fig. |l5|, with all parameters un- 
changed except huiph — 90meV. At 10K and 20K, a split- 
ting of the principal dos peaks is just resolvable. This 
splitting is due to the presence of combined CDW and 
superconducting order, but it does not simply mean that 
the two phases have separate spectral peaks. Thus, at 
(n, 0), there is a single peak, at the position given by 
Eq. m. Instead, the two gaps have different dispersion, 
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with the superconducting gap always at 6f and the CDW 
gap dispersing away from tp (see feature C in Fig. ||b). 
Thus, at (it, 0) the gap is the (Euclidean) sum of the 
CDW and superconducting gap, while near (tt/2, tt/2) it 
is only the latter. The splitting of the peaks in Fig. Ore- 
fleets this anisotropy. This splitting is essentially lost by 
30K, due to thermal broadening, and would not survive 
much disorder broadening. Note that the gap fills in with 
increasing temperature, without significantly closing. In- 
deed, since superconductivity disappears at a lower tem- 
perature, the gap actually appears to increase at higher 
temperatures. The low-temperature splitting increases 
with increased underdoping (Fig. |l6| ), and hence might 
ultimately be observable experimentally. However, the 
data of Fig. [l4| correspond to T c = 69K, somewhat lower 
than any reported tunneling data. Moreover, the calcu- 
lations are for an s-wave superconductor, and the smaller 
gap could be very different, or even absent, for a d-wave 
superconductor. _ 

Thus, the model agrees with experiments in finding 
gap-like features above T c , but does not reproduce the 
rearrangement of features below T c into a sharp peak, 
dip, plus hump. A suggested explanation for these fea- 
tures is given in the following Section. 

Figure [l6] shows the evolution of the low-temperature 
tunneling spectra with doping, in the underdoped regime. 
At optimal doping (lowest curve) the gap is symmetrical 
and pure superconducting. As the material is succes- 
sively underdoped, the peaks split, with the pseudogap 
peak growing and the superconducting peak shrinking, 
scaling with T c , down to x ~ 0.06 where superconductiv- 
ity disappears. 



transfer gap of ~ 1.6eV, and a residual dispersion sugges- 
tive of the flux phase. Figure [l?] shows that both these 
features are qualitatively present in the model. From 
Eq. ^l], it can be seen that near half filling, when t' — > 0, 
the CDW gap opens up over the full Fermi surface, while 
the excess G\ gap ensures that the largest gap is present 
near the VHS. The chief differences from experiment are 
(a) the gap evolves discontinuously from Go to Go + G\ 
rather than smoothly, (b) the overall bandwidth is too 
large, since the model has not taken into account that 
correlations renormalize the bandwidth from ~ 8t to 
~ 2JE3, and (c) the model cannot reproduce the mag- 
nitude of the (electronic) Mott gap, ~ 1.6eV. Indeed, 
the full gap is limited to < 2hu> P h- 
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FIG. 17. Band dispersion near half filling Bi-2212 at 10K, 
using parameters of FigM 
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FIG. 16. Tunneling dos for underdoped Bi-2212 at 10K, as 
a function of doping, using parameters of Fig.|§| From top 
to bottom, the doping is x = 0.004, 0.04, 0.08, 0.124, 0.169, 
0.218, 0.274, and x c = 0.34. All curves are offset for clarity 
(all approximately coincide at eV = 200meV). 

The lowest doped state, close to half filling, is very in- 
teresting. For this state the present model is expected 
to be least accurate. The state is dominated by an elec- 
tronic instability, the Mott transition. There is a charge- 



E. Peak, Dip, and Hump 

The photoemission in Bi-2212 near (tt, 0) has a charac- 
teristic shape in the superconducting state which is diffi- 
cult to reproduce in the current model. In the pseudogap 
state above T c , there is a very broad hump, consistent 
with a CDW gap with large broadening T, as in Section 
II. D. In the superconducting state, the hump remains, 
but the low-energy side of the hump is transformed into 
a very sharp feature at energy Ai, with a dip at energy 
~ 2Ai. In the CDW superconductor, there is not ex- 
pected to be such a two-peaked structure along (n, 0), 
but a single peak given by Eq. |l|. However, as noted in 
Section II. D, the broadening can have a strong energy 
dependence, with a long lifetime in the superconducting 
state cut off near 2Acj. In the absence of a full the- 
ory, it is not clear which gap A represents, but since this 
is specifically a coherence effect, it seems reasonable to 
assume that it should be approximately the supercon- 
ducting Afc. Since A& < A t , one can get a behavior 
at low temperatures which is closer to experiment. Fig- 
ure [l8] illustrates the three possible behaviors, depending 
on whether 2A& is <<, ~<, or > A t : one can get hump 
plus dip (dotted line), hump plus peak (dashed line), or 
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peak plus dip (solid line) . In the last case, the dip is not 
necessarily at twice the peak in energy. 
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FIG. 18. (it, 0) photoemission spectra for Bi-2212 at 
T=1K, using parameters of Fig.| (x=0.22), except Ag = 0.446 
(dotted line), 0.385 (dashed line), and 0.325eV (solid line). 
Broadening given by Eq. E3l 
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FIG. 19. (-7T, 0) photoemission spectra for Bi-2212, using 
parameters of Figji (x=0.22), except inhomogeneously broad- 
ened (averaged over 20 Ag values). Solid line: T = IK, dashed 
line: T = 70K. 

It seems clear, however, that the hump is associated 
with the pseudogap, and the peak with the supercon- 
ducting gap. We suggest that the experimental result 
may represent an effect of phase separation, which is not 
well captured by our model. To test this hypothesis, we 
developed a simple model, assuming that the pseudogap 
feature is inhomogeneously broadened. Specifically, we 
assume that Ag is spatially inhomogeneous, and that re- 
gions of different Ag evolve independently. For each re- 
gion, lowering the temperature will produce a spectrum 
like one of those in Fig. [l^. In the regions of large Ag, 
superconductivity has little effect, leaving the hump un- 
changed. However, when Ag is small, at low T there will 
be a sharp superconducting peak plus well-defined dip. 



Figure 19 shows the resulting spectra, both just above 
the onset of superconductivity, at T=70K (dashed line), 
and at low temperature, T=1K (solid line). Each spec- 
trum is the superposition of 20 individual spectra, each 
with 



r = 



r x , if \E\ < 2A fe 
T2, otherwise, 



(23) 



with Ti — 2meV, T2 = lOmeV. The general features 
of the experiment are clearly reproduced, although the 
sharp peak is not resolution limited, due to the inhomo- 
geneous broadening. 

This interpretation of the sharp peak as a measure of 
Afc leads to a new problem, since the experiments imply 
that this Afe increases with decreasing x. To reproduce 
this behavior in the phase diagram would require that Aa 
varies with doping - the same effect needed to explain the 
LSCO phase diagram, Fig. [| Such a doping dependence 
would be consistent with magnetic fluctuation-induced 
superconductivity. However, it must be kept in mind 
that an alternative interpretation, such as the solid line 
in Fig. [ill was ruled out because the dip feature falls 
at esactly 2A X . There exist some data (e.g., Fig. 1 of 
Ref.Ej) in which the dip falls at an energy substantially 
less than 2A X . This discrepancy must be clarified before 
a definitive model for the sharp peak can be established. 



F. The Incredible Shrinking Fermi Surface 

Norman, et al.El have shown that the Fermi surface 
in the pseudogap phase has a remarkable temperature 
dependence: there is a full, large Fermi surface above 
the pseudogap transition, T p , but as T is reduced be- 
low T p , the Fermi surface gradually collapses - vanishing 
first near (n, 0) and then over a larger angular range be- 
fore ultimately being reduced to a point along the line 
c x = c y . This low temperature limit is consistent with 
a d-wave superconductor, but the intermediate tempera- 
ture regime is not: the gap should open at T c everywhere 
on the Fermi surface, with the magnitude of the gap pro- 
portional to (c x — c v ). Fluctuation broadening could pro- 
duce this effect by limiting the pseudogap to those parts 
of the FerraL surface where it is larger than the energy 
uncertaintyE^TTJ. However, it can also be interpreted in 
terms of a VH nesting gap. 

Figure 20 shows the evolution of the Fermi surface with 
Go, for A = G\ = 0. In a pure CDW state, the Fermi 
surface can be calculated analytically, Appendix C. For 



finite Go, the surface is an ellipse, with the inner half de- 
rived from the original Fermi surface and the outer half 
a ghost Fermi surface, zone folded by Q. The latter has 
considerably lower intensity (inset), and will be further 
reduced by fluctuation effects and stripe effects; it is not 
seen experimentally. On the other hand, the truncated 
original surface has a large coherence factor (> 0.5), and 
bears a strong resemblance to the experimental data. Be- 
low T c , the Gj's shrink as A grows, so the opening of the 
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superconducting gap takes over. Thus, for a d-wave su- 
perconductor, the residual points of zero gap would lie 
along the original Fermi surface, as observed, and not at 
(7r/2,7r/2). When G x ^ 0, Fig. |l|, the resemblance to 
experiment is even closer: the Fermi surface shrinks down 
with essentially no change of shape, from the pseudogap 
onset to the superconducting T c . 
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FIG. 20. Fermi surface in CDW phase, for r = -0.38, Go 
= 20, 40, 80, 160, and 180 meV, for progressively smaller 
ellipses. Dotted segments = ghost Fermi surfaces, with co- 
herence factors < 0.5. Inset: coherence factors along the 
diagonal (k x — k y ). 
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locus found in the pseudogap phase. These problems will 
be discussed further in the following section. 



G. Minimum Gap Locus 

In order to recover a 'Fermi surface' in the presence of 
a (pseudo)gap, Ding, et al.E3 have introduced the con- 
cept of a 'locus of minimum gap'. For a series of cuts 
in the Brillouin zone, they define the minimum gap - 
the fc-point at which the photoemission gap most closely 
approaches the Fermi level. The gap is defined as the 
leading edge of the photoemission pattern, but the com- 
parison with tunnelinga suggests that the peak position 
should work as well. The Fermi surface would then cor- 
respond to a locus of zero-gap excitations. In Appendix 
C, analytical expressions for the Fermi surface and mini- 
mum gap locus (actually, the maximum in the dispersion 
of the lower band) are derived, which are plotted as lines 
in Fig. g2[ The solid lines are the Fermi surface, the 
dashed lines the minimum gap loci, and the dot-dashed 
and dotted lines their ghostly counterparts. Neglecting 
the ghost surfaces, the-xesulting dispersion resembles the 
experimental findinga£3. Note in particular that the locus 
of minimum gap does not pass through the (7 r, 0) point, 
even though is pinned at the VHC (see Eq. C4). This 
can readily be understood: since the gap is largest at 
(77, 0), the scans find a smaller gap at some k x away from 
this point. Knowing the full dispersion E(k), a much 
fuller test of the theory can in principle be carried out. 



(O.tt) 



(7T, 7t) 
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FIG. 21. Fermi surfaces in CDW phase, for self-consistent 
parameters. Lower (upper) set of lines: r = -0.4 (-0.16), Xg 
— 0.50 (0.25) eV. Dotted lines = in normal state above T v \ 
solid lines = in pseudogap state, just above T c ; dashed lines = 
intermediate temperatures. Only segments of Fermi surface 
with coherence factors > 0.1 are shown. 



There is one possibly significant difference: as the nor- 
mal state above T p is approached, Fermi surfaces should 
approach the VHS's at (it, 0), which is not clearly seen 
in experiment. A related problem is the minimum gap 
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FIG. 22. Loci of minimum gap (dashed and dotted lines) 
and zero gap (solid and dot-dashed lines). 

While the problem of the minimum gap locus can be 
satisfactorily explained in the pinned VHS model, the 
problem remains that in the normal state above the pseu- 
dogap, the Fermi surface found in photoemission does not 
seem to intersect (7r, 0). It must be recalled that the VHS 
is the ultimate hot spot, so local residual disorder is likely 
to broaden or split the quasiparticle spectrum near the 
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VHS even at temperatures considerably in excess of the 
pseudogap temperature. An indication fpt such smearing 
is the disagreement between the ArgonneE3 and StanfordEJ 
groups on the shape of the Fermi surface near the VHS 
in underdoped cuprates. Hence we prefer to infer the 
position of the VHC from data in the superconducting 
state, when line broadening effects are smallest. If the 
VHC is within 5meV of in the superconducting state 
(Fig. |l^), it should be even closer in the normal state, 
since the superconducting transition tends to shift the 
chemical potential away from the VHSE^I. 

H. Spectral Weight Shift 

Shen, et alB have recently found significant shifts of 
spectral weight on cooling the sample from the pseudo- 
gap phase above T c to the superconducting phase below. 
These are exactly the sort of changes one would expect in 
the present model, since the CDW and superconductivity 
have very different gap structures, and the CDW gap is 
significantly reduced below T c , Fig. ^. However, experi- 
mentally these changes are found to occur over an energy 
range extending to ~ SOOmeV away from ep. This ex- 
tended range is presumably telling us about the relevant 
pairing bosons, but since these energies are greater than 
Tiujph, strong coupling calculations are required, which 
are beyond the scope of the present paper. 

IV. CONCLUSIONS 

We have taken a simple pinned Van Hove Ansatz for 
the striped pseudogap phase in the cuprates, and ana- 
lyzed the predicted tunneling and photocmission spec- 
tra, comparing them with experimental data. Highlitcs 
of our results include: (1) We explain the experimental 
observation that the tunneling peaks coincide with the 
(jr, 0) photoemission dispersion, and show that careful 
measurements of this effect can provide direct evidence 
of pinning of the Fermi level to the VHC over an extended 
doping range. (2) In turn, the fact that tunneling shows 
a well-defined gap-like feature confirms that the (ir, 0) 
dispersion also has two branches - that is, that the_pseu- 
dogap is associated with some form of VHS nestingtd. (3) 
The tunneling gap has a characteristic asymmetry which 
vanishes at optimal doping; this is evidence that opti- 
mal doping is that point-at which the Fermi level exactly 
coincides with the VHSE3. 

(4) By plotting the doping dependence of the photoe- 
mission hump feature, the resulting pseudogap phase dia- 
gram for Bi-2212 is in good agreement with similar phase 
diagrams for LSCO and YBCO, derived from transport 
measurements. (5) The rearrangement of spectral weight 
seen below T c can be interpreted as a generation of sep- 
arate pseudogap (hump feature) and superconducting 
(sharp feature) peaks with the dip between them due 



to superconducting coherence effects. However, this re- 
quires inhomogeneity not included in the simple Ansatz. 
(6) Finally, a number of specific features observed in 
the photoemission receive a natural explanation in this 
model, including the shrinking Fermi surface and the lo- 
cus of minimum gap. 

A shortcoming of the model is the prediction of a split- 
ting of the tunneling gap, which has not yet been ob- 
served. However, we note that (a) this may be due to 
the assumption of s-wave superconductivity, and (b) even 
the assumption of a pure d-wave superconductor, with no 
pseudogap, would lead to a split peak in tunneling unless 
the Fermi level were pinned to the VHS. 

V. DISCUSSION 
A. Alternative Pseudogap Scenarios 

The present interpretation of the pseudogap in terms of 
Van Hove n es t in g is in agreement with some early (1990) 
calculations! 6 jC3 which predated any of the experimen- 
tal observations of the 'spin gap' or pseudpgap. How- 
ever, ona.of these theories involves an SDWEi the other 
a CDWEj, while more recent calculationsE3ta involve a 
flux phase near half filling. Hence, the essential feature is 
the Van Hove nesting; the particular instability depends 
on the details of the electron-electron interaction. 

It is convenient to classify these instabilities in terms of 
the SO (6) scenario. As discussed in Section II. B, the in- 
stabilities fall into two broad classes: nesting and pairing. 
We expect that the phase diagram, Fig. [j], of the compe- 
tition between pseudogap and superconductivity, would 
not be greatly changed if the pseudogap represents any 
nesting instability, as long as the physics is dominated 
by proximity to a VHS. This would hold for any Hub- 
bard or tJ model, in which the VHS is at half filling, 
and more particularly for any model which introduces a 
higher order hopping parameter (t 1 or too) to move the 
VHS toward finite positive x (hole doping) - i.e., closer 
to optimal doping. 

Among the nesting theories of the pseudogap, jne in- 
clude Schrieffer's spin bagE3 and Zhang's J=jO(5)El, for 
which the_instability is an SDW, Laughlir£3's and Wen 
and LeeEfs flux phase models, and Klemm's CDW 
modelEJ. Within an SO (6) symmetric model, all would 
produce the same pseudogap phase diagram. Hence, 
the question of which one actually produces the ground 
state depends on which symmetry-breaking operators are 
present in the physical cuprates. 

Following the discussion in Section II. B, it would be 
more difficult to reproduce the pseuodgap phase dia- 
gram in terms of a pairing instability - i.e., assuming that 
the pseudogap is a signature of superconducting fluctua- 
tions, a precursor of real space pairing. We believe that 
this class of models will have difficulty explaining striped 
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phases, and the seemingly smooth extrapolation of the 
pseudogap to the SCOC dispersion at half filling. 



B. Distinguishing the Nesting Instabilities 

There is a natural generalization!! of Eq. [I] to SO (6): 



\ 
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(24) 



where the sum is over all twelve instabilities of both su- 
perspins. Thus, at (n, 0) the pseudogap depends only on 
the vector sum of the individual gaps. For t' = 0, Eq. |24| 
holds over the full Fermi surface. Hence, the instabilities 
can only be distinguished by their dispersion away from 
(tt.O). 

In attempting to distinguish different models, it must 
be kept in mind that this secondary dispersion is likely to 
be model dependent. For example, d-wave superconduc- 
tivity would have different dispersion depending whether 
it was phonon-induced or an electronic instability. In 
particular, the dispersion will be sensitive to the nature 
of the boson mediating superconductivity, in particular 
to the cutoff, lj c . 

A second example is instructive. In the antiferromag- 
netic state at half filling there is a competition between 
the Neel phase (commensurate SDW) and the flux phase. 
In a Hartree-Fock calculation, both phases give rise to 
similar gaps at (7r,0), but very different dispersion in 
other directions. In particular, the Neel gap is approxi- 
mately isotropic, while the flux-^phase has zero gap near 
(tt/2, it/2). (See Fig. 32 of RefEO and the associated dis- 
cussion for references to earlier literature.) Having the 
larger average gap, the Neel phase is-, more stable, but 
the experimental dispersion of SCOQia, and it was spec- 
ulated that the Neel phase dispersion is modified by fluc- 
tuation or correlation effects. This seems to be the case, 
since calculations of the dispersion of one hole in an anti- 
ferromagnet are capable of reproducing the experimental 
dispersionEj. Thus, while the theoretical interpretation 
has proven quite involved, the dispersion remains sim- 
ple: it has the preiodicity of the Neel superlattice, with 
maximal gap associated with VHS splitting. 

Despite these caveats, we are exploring an SO(6)- 
generalized version of the pBF model, to test how sensi- 
tive the phase diagram and tunneling and photoemission 
spectra are to the particular forms of nesting or pairing 
instabilities. 



C. Doping Dependence of Aa 

There were a few tantalizing hints that the supercon- 
ducting coupling parameter Aa is doping dependent: in- 
creasing as x is reduced below optimal doping. This was 
needed to explain the LSCO phase diagram, Fig. ||, and 



the doping dependence of the sharp photoemission peak 
below T c , Section III.E. This could have a simple expla- 
nation - e.g., the bandwidth is reduced by correlation 
effects from ~ At (half bandwidth) at optimal doping to 
~ 2 J near half filling, and this reduction is neglected in 
the pBF model. Alternatively, it could signal the impor- 
tance of magnetic fluctuations. ._. 

In the three-band slave boson model calculations, the 
stripes were interpreted as a crossover from a flux phase 
near half filling to a CDW phase near optimal doping. In 
SO (6), this is a crossover between the two different super- 
spin multiplets, and is accompanied by a crossover from 
d-wave to s-wave superconductivity. This would be con- 
sistent with hints of a more s- wave-like superconducting 
transition in overdoped materials. 

We would like to thank NATO for enabling A.M. 
Gabovich to visit us and discuss his work. Publication 
744 of the Barnett Institute. 



APPENDIX A: PHASE DIAGRAM FOR LSCO 
AND YBCO 

We briefly discuss the procedure for determining the 
phase diagram of Figs. [l||3| and ||, to indicate the role of 
the various parameters. 

1. Choice of r. Using Eq. S, a value r = 2t' /t = 
—0.38 is necessary for the VHS to fall at optimal dop- 
ing x c — 0J6 in LSCO. On the other hand, there is a 
suggestionEJ that r in LSCO may be as large as -0.6, 
which corresponds to x c = 0.274. 

In YBCO, band structure calculations suggest r = 
—0.9 - remarkably close to the critical value -1.0 for 
one-dimensional behavior. This in turn corresponds to a 
doping x c = 0.54. Note however, that the band structure 
calculations also find a significant interplane coupling be- 
tween the two Cu02 planes in a unit cell, such that the 
symmetrical combination is essentially undoped. In this 
case, for the antisymmetrical combination Fermi surface 
to be at the VHS would require an average doping per 
layer of 0.54/2=0.27 holes. This is consistent with the op- 
timal doping estimated from the Uemura plotEJ, which 
finds that x c /m is about twice as large for YBCO as for 
LSCO. Assuming that both materials have the same ef- 
fective mass m, this would give a broad peak in x c near 
0.32 for YBCO. 

The situation for Bi-2212 is slightly more complicated. 
The Uemura plot suggests an x c similar to that ituYBCO, 
while photoemission is consistent with r ~ — 0.9lj. How- 
ever, there is no evidence in photoemission for the inter- 
layer splitting of the Fermi surface. We will postulate 
that a splitting similar to that found in YBCO is locally 
present, but is obscured by a very short c-axis correlation 
length. 

2. Optimal Doping. In YBCO, the pseudogap ap- 
pears to vanish very close to optimal doping, as soon as 
the superconducting T c exceeds the pseudogap T p . This 
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is a natural consequence of the BF model, and so is a con- 
venient point to fix for the phase diagram. We proceed as 
follows. For a given r, we set the Fermi level at the VHS 
and calculate the electron-phonon coupling energies, Aa 
and A G which give A k = G k = 0+ at T = 90K. With 
this choice of A's, T c — T p = 90K at this doping, and 
at any larger doping the pseudogap will be suppressed to 
zero by T c . 
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FIG. 23. Allowed values of A vs u> p h to fix T, 



at optimal doping. Solid (dashed) line = 
(-0.9); dotdashed (dotted) line = Aa for t 



p c = 90K 
A<3 for t — -0.6 
= -0.6 (-0.9). 



The only free parameter is uj p h- In Fig. [2^, we plot the 
resulting A's vs Lu p h for two choices of r. As expected, for 
a larger w p h, a smaller A is needed, but for io p h > 30meV, 
the value of A varies only weakly. The dependence of the 
A's on t is what would be expected from Fig. |f nesting 
is strongest for r = 0, so the necessary Ag increases with 
increasing r, while pairing becomes easier as the bands 
become more one-dimensional, r — > 1. 




FIG. 24. Gap ratio vs ui p h using values of A determined 
from Fig. Solid (dashed) line = ratios for superconducting 
gap for r = -0.6 (-0.9); dotdashed (dotted) line = ratios for 
CDW gap for r = -0.6 (-0.9). 



in particular: the gap ratio, the magnitude of the limit- 
ing pseudogap at zero doping, and the superconducting 
onset. 

In weak coupling BCS theory, the gap ratio 2A(T = 
0)/fcsT c = 3.53, and increases for stronger coupling. Fig- 
ure [24] shows the variation of this ratio with uj p h , where 
the A's are fixed by Fig. |2^, as discussed above. The ra- 
tio for the superconducting phase is somewhat enhanced 
from the BCS value, and depends only weakly on uj p h- 
The ratio can be much larger and more w p ^-dependent 
for the pseudogap. Note that in Fig. G], where Aa = Ag is 
assumed, the two gap ratios are nearly equal. (The prob- 
lem with Fig. |l| is that the,-critical doping is too low.) 
Experimentally, Oda, £t alia find a value 4-5, whereas 
the data of Ding, et al.Eil yield ~ 8. Thus, the model can 
account for the former value, but not the latter. If the 
correct ratio turns out to be ~ 8,,-this may be evidence 
for pairbreaking effects limiting TJI3. 
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FIG. 25. CDW onset T v at zero doping vs u p h using values 
of A determined from Fig. |23| Dotdashed (dotted) line = T p 
for CDW gap for r = -0.6 (-0.9). 

Figure ^ shows how the limiting pseudogap T p {r = 0) 
depends on uj p h (in view of the uncertainty in the gap ra- 
tio, we have chosen to work with the critical temperature 
rather than the gap). 
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3. Choosing uj p t If we assume t — 0.25eV is known, 
and that Aa and Ag are doping independent, there is 
only one remaining free parameter, u p h, to adjust to fit 
the rest of the phase diagram. We explore three features 



FIG. 26. Onset of superconductivity x vs ui p h using values 
of A determined from Fig. |2^. Dotdashed (dotted) line = x 
for CDW gap for r = -0.6 (-0.9). 
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Another constraint is the superconducting onset x m i n , 
the doping at which superconductivity first appears. 
This is - 0.05 in LSCO; for YBCO and Bi-2212, it is 
hard to specify the precise planar hole doping. Adjusting 
Xmin is somewhat in conflict with the other constraints, 
since a large difference Ac ~ Aa leads to a larger x m i n , 

Fig. El. 

Given these conflicting constraints, the parameters of 
Fig. H seem to provide the best overall solution. 



APPENDIX B: ANALYSIS OF D-WAVE GAPS 



using Sy = 1 — c 2 and g — G^/2t, in addition to the 
notation introduced in Appendix B. The Fermi surface 
shrinks to the point (ir/2, ir/2) when g = \t\. 

The locus of minimal gap can be found as in Appendix 
B. Its exact position depends on the orientation of the 
scan. Here, we fix k y and scan k x . The top of the E^- 
band is then at 



Z X Cy i 



In particular, for c y = 1 



\Tgc y \ 



1 - T 2 C 2 



(C3) 



For a purely superconducting gap (Go = G\ = 0) the 
eigenvalues of Eq. ^ simplify: 
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) 2 + A|. 
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For a d-wave superconductor, = Qj:Aq(c x — c y ). The 
tunneling peak corresponds to the maximum of (or 
minimum of E + ). From Fig. [H^, it is sufficient to evaluate 
dE-ldc x when c y — ±1 (the +1 solution corresponds to 
overdoping, -1 to underdoping). In terms of the scaled 
variables r = 2t'/t, S — Ao/2£, and fi = ep/2t, the 
maximum occurs at 



± c x = 



-(l±r)(l±ft) 
5 2 + {1±t) 2 ' 

and is E max = Apu with 

,E max ,2 6 2 (r + fi±2) 2 



2t 



5 2 + (1±t) 2 ' 
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which give the dashed lines in Fig. [H]. The photoemission 
gap Ape (solid line in Fig. [ll]) is simply the value of £L 
at (tt,0): 



(^) 2 = (r- A ) 2 +^. 



(B4) 



When the Fermi level coincides with the VHS, p, = t, 
and Ape = Atu = 2Aq, to lowest order in S 2 . 



APPENDIX C: PROPERTIES OF THE CDW 

When superconductivity is absent, the Fermi surface 
can be found from solutions to 



r 2 



(CI) 



where eg = — ep- When the Fermi level is pinned at 
the VHS, tp = At' , this reduces to 



-(l-r 2 )c y ± ^r 2 Sy-g 2 {\-T 2 c 2 y) 

l - T 2 C 2 y : 



(C2) 



-1 



\rg\ 



(C4) 



Note that, even though the Fermi level is pinned at the 
VHC, the locus of minimum gap does not pass through 
(tt.O). 

* : On leave of absence from Inst, of Atomic Physics, 
Bucharest, Romania 
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